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ABSTRACT 

The Stephan's Quintet (hereafter SQ) is a template source to study the impact of galaxies interaction on the 
physical state and energetics of their gas. We report on IRAM single-dish CO observations of the SQ compact 
group of galaxies. These observations follow up the Spitzer discovery of bright mid-IR H2 rotational line 
emission (L(H2) ~ lO-'^ W) from warm (10^""' K) molecular gas, associated with a 30 kpc long shock between 
a galaxy, NGC 7318b, and NGC 7319's tidal arm. We detect CO(l-O), (2-1) and (3-2) Une emission in the 
inter-galactic medium (IGM) with complex profiles, spanning a velocity range of sa 1000 km s"'. The spectra 
exhibit the pre-shock recession velocities of the two colliding gas systems (5700 and 6700 km s"'), but also 
intermediate velocities. This shows that much of the molecular gas has formed out of diffuse gas accelerated 
by the galaxy-tidal arm collision. CO emission is also detected in a bridge feature that connects the shock to 
the Seyfert member of the group, NGC 7319, and in the northern star forming region, SQ-A, where a new 
velocity component is identified at 6900 kms"', in addition to the two velocity components already known. 
Assuming a Galactic CO(l-O) emission to H2 mass conversion factor, a total H2 mass of « 5 x 10^ M© is 
detected in the shock. The ratio between the warm H2 mass derived from Spitzer spectroscopy, and the H2 
mass derived from CO fluxes is « 0.3 in the IGM of SQ, which is 10-100 times higher than in star-forming 
galaxies. The molecular gas carries a large fraction of the gas kinetic energy involved in the collision, meaning 
that this energy has not been thermalized yet. The kinetic energy of the H2 gas derived from CO observations 
is comparable to that of the warm H2 gas from Spitzer spectroscopy, and a factor « 5 greater than the thermal 
energy of the hot plasma heated by the collision. In the shock and bridge regions, the ratio of the PAH-to-CO 
surface luminosities, commonly used to measure the star formation efficiency of the H2 gas, is lower (up to a 
factor 75) than the observed values in star-forming galaxies. We suggest that turbulence fed by the galaxy-tidal 
arm collision maintains a high heating rate within the H2 gas. This interpretation implies that the velocity 
dispersion on the scale of giant molecular clouds in SQ is one order of magnitude larger than the Galactic 
value. The high amplitude of turbulence may explain why this gas is not forming stars efficiently. 
Subject headings: Galaxies: clusters: individual: Stephan's Quintet - galaxies: interactions - galaxies: ISM - 
intergalactic medium 



1. INTRODUCTION 

Galaxy interactions represent an important agent of galaxy 
evolution, involving the injection of huge amounts of kinetic 
energy in the interstellar medium. Galaxy collisions or in- 
teractions of galaxies with the inter-galactic medium (hence- 
forth IGM) may trigger bursts of star formation detected in 
the infra-red (IR) (e.g. Joseph & Wright 1985; Harwit et al. 
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1987; Kim et al. 1995; Xu et al. 1999). Star formation pro- 
ceeds where and when the kinetic energy of the gas is being 
radiated, allowing gas to cool and condense. To make head- 
way in our understanding of galaxy evolution and, in partic- 
ular, the star formation history of the universe, it is crucial to 
understand how this kinetic energy is dissipated, and how this 
dissipation affects the passage from molecular gas to stars. 

Stephan's Quintet (Hickson Compact Group HCG 92), 
hereafter SQ, is an ideal laboratory to study these physical 
processes. SQ is a group of five interacting galaxies, with a 
complex dynamical history (see Xu 2006, for a review and 
references therein). The most striking feature of the group 
is that a galaxy-scale (« 15 x 35 kpc") shock is created by 
an intruding galaxy, NGC 7318b which is colliding into the 
intra-group gas at a relative velocity of « 1 000 km s"' (see the 
background image in Figure 1). This IGM material is likely to 
be gas that has been tidally stripped from NGC 73 19 by a past 
interaction with another galaxy, NGC 7320c (see Moles et al. 
1997; Sulentic et al. 2001, for a plausible dynamical sce- 
nario). These ideas have been further explored recently with 
numerical simulations. Both N-body (Renaud et al. 2010) and 
hydrodynamical models (Hwang et al. 2011) of the Quintet 
suggest a sequence of serial collisions between various group 
members culminating in the current high-speed shock-heating 
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of gas as NGC 7318b strikes the pre-existing tidal debris from 
NGC7319. 

A ridge of X-ray (Trinchieri et al. 2003, 2005) and strong 
radio synchrotron (e.g., Sulentic et al. 2001; Williams et al. 
2002) emission from the hot (5-7 x lO^K) post-shock plasma 
is associated with the group-wide shock. Excitation diagnos- 
tics from optical (Xu et al. 2003) and mid-lR (Cluver et al. 
2010) emission lines also confirm the presence of shocked 
gas in the IGM. Spitzer mid-lR observations reveal that 
the mid-lR spectrum in the shock structure is dominated 
by the rotational line emission of molecular hydrogen, H2 
( Apple ton et al. 2006; Cluver etal. 2010). This denotes the 
presence of large amounts (w K^Mq) of warm (10^- lO-'K) 
H2, co-spatial with the X-ray emitting plasma of the group- 
wide shock. The H2 emission is extended (see the H2 red 
contours on Fig 1) along the ridge feature, towards the Seyfert 
galaxy NGC 7319 (in the so-called bridge), as well as in the 
northern intergalactic starburst SQ-A (Xu et al. 1999). Over 
the main shock area observed by Spitzer, the H2 luminos- 
ity (L(H2) w 2.6 X 10*^ Lq, summed over the O-OS(O) to S(5) 
lines) is larger by a factor of three than that of the X-ray 
emission (integrated over 0.001 - lOkeV) from the same re- 
gion. With the high resolution {TZ = 600) of Spitzer IRS, the 
17 pm H2 S(l) hne is resolved with a FWHM of 870 km s"' 
(Appleton et al. 2006). 

Guillard et al. (2009) proposed a model for the formation of 
molecular gas in the ridge, considering the collision of multi- 
phase gas. In their scenario, H2 is formed from the cooling of 
shocked atomic gas. The shock velocity depends on the pre- 
shock gas density: it is high (« 600- 800km s~') in the tenu- 
ous («H < 0.02 cm"-') gas, producing the X-ray emission in the 
ridge, and much lower in the dense gas. Gas at pre-shock den- 
sities «H > 0.2 cm~^ experience shock velocities small enough 
(Vs < 200 km s"') to retain most of its dust, cool, and become 
molecular within a few million years. We proposed that the 
H2 emission is associated with the dissipation of the kinetic 
energy of the galaxy collision at sub-parsec scales, through 
supersonic turbulence in the molecular gas. The mid-lR spec- 
tral energy distribution (SED) of the H2 emission can be ex- 
plained if the dense («h > lO^cm"^) molecular gas is continu- 
ously processed by low-velocity (5-20 km s"') MHD shocks. 

Since the mid-lR rotational H2 line emission only samples 
warm molecular gas with temperatures « 10^- lO"* K, our cen- 
sus of the molecular gas in the SQ shock is still incomplete. 
The discovery of powerful mid-lR H2 line emission in SQ led 
us to search for CO line emission associated with the warm 
H2 in the shock. In galaxies, the warm H2 is a small fraction 
(0.4-4 %) of the total mass of molecular gas (Roussel et al. 
2007), which is too cold (< 100 K) to be seen in mid-lR 
H2 emission, and usually traced with CO. This paper reports 
on the detection of CO line emission in the IGM of the SQ 
group, and investigates the physical and dynamical state of 
the molecular gas in SQ, as well as the connection between 
shock and star formation in the group. 

After reviewing the past CO observations of SQ (Sect. 2), 
we present technical details about our observations with the 
IRAM 30-meter and the APEX telescopes (Sect. 3). Then we 
describe the kinematics (Sect. 4), the distribution and mass 
(Sect. 5), the excitation and energetics (Sect. 6) of the molec- 
ular gas in the IGM of the group. Sect. 7 discusses the origin 
of the CO gas, its kinematics, and the efficiency of star for- 
mation in the group. Conclusions and final remarks are given 



in Sect. 8. 

Throughout this paper we assume the distance to the 
SQ group to be 94 Mpc (with a Hubble constant of 70 
km s~' Mpc~~) and a systemic velocity for the group as a 
whole of 6600 km s~'. At this distance, 10 arcsec = 4.56kpc. 

2. PREVIOUS CO OBSERVATIONS IN THE SQ GROUP 

Before the discovery of warm H2 in the IGM of the SQ 
group, previous observations were mostly concentrated on 
the Seyfert galaxy NGC 7319, on the southern tidal fea- 
tures, or on the northern starburst region (SQ-A). The first 
millimetre CO observations in SQ detected cold molecular 
gas in NGC 7319 (Yun et al. 1997; Verdes-Montenegro et al. 
1998; Leon et al. 1998), which contains 4.8 x lO'^ M© 
Smith & Struck (2001). Using the BIMA" interferometer, 
Gao & Xu (2000); Petitpas & Taylor (2005) showed that 3/4 
of the CO gas in NGC 7319 is lying outside its nucleus, 
mostly in an extended northern region (see green contours on 
Figure 1). 

Molecular gas outside galactic disks, associated with 
the IGM starburst SQ-A, was first detected by Gao & Xu 

(2000) using BIMA. Then, single dish observations by 
Smith & Struck (2001) and Lisenfeld et al. (2002) confirmed 
this result. Two velocity components were detected, one cen- 
tred at 6000 km s"', and the other at 6700 km s"'. 
These velocities match the redshifts of the two Hi gas 
systems found in the same region (Sulentic et al. 2001; 
Williams et al. 2002), which correspond to the intruder 
galaxy, NGC 7318b, and NGC 7319's tidal tail in the IGM, 
respectively. Lisenfeld et al. (2002) found that there is more 
molecular gas (3.1 x lO'' M©) in SQ-A than Hi gas (1.6 x 
lO"* Mq). 

Past observations with the BIMA interferometer did not 
lead to a CO detection in the shock itself (Gao & Xu 2000). 
The single-dish observations by Lisenfeld et al. (2002) using 
the IRAM 30m telescope partially overlap the shock struc- 
ture and show a 3(t detection of CO in the northern region 
of the X-ray ridge, but their observations do not probe the 
entire shock region. We note that they include this emission 
in their SQ-A CO flux. Gao & Xu (2000); Smith & Struck 

(2001) ; Petitpas & Taylor (2005) reported CO emission close 
to NGC 7318b nucleus and from several other regions in the 
group. 

Another interesting feature of SQ is the presence of a 
co-rich (Braine et al. 2001; Lisenfeld et al. 2002, 2004, 7 x 
10** Mq) region, SQ-B, located within the young Hi tidal tail 
(see Sulentic et al. 2001), to the south of NGC 7319 (outside 
the field of view of Figure 1). This region is also bright in the 
mid-lR (Xu et al. 1999; Guillard et al. 2010), Ha (Arp 1973) 
and UV (Xuetal. 2005). The H2-to-Hl gas mass ratio of 
0.5 indicates that SQ-B is a tidal dwarf galaxy (TDG) candi- 
date, i.e. a small galaxy which is in the process of formation 
from the material of the NGC 7319's tidal tail (Boquien et al. 
2009). The gas metallicity in SQ-A and SQ-B is slightly 
higher than solar (Xu et al. 2003; Lisenfeld et al. 2004). This 
is a strong indication that the gas in these regions has been 
pulled out from the inner part of a galaxy disk (or maybe sev- 
eral disks) by tidal interactions. TDGs will not be further dis- 
cussed, as we mostly concentrate on the shock, the bridge and 
SQ-A in this paper 

" Berkeley Illinois Maryland Association, 
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Table 1 

CO observations log 



Position 


RA (J2000) 


DEC (J2000) 


integration time [min] 
CO(l-O) CO(2-l) CO(3-2) 


SQ-A" 


22:35:58.88 


+33:58:50.69 


122 


58 


66 


Ridge 1 


22:35:59.85 


+33:58:16.55 


311 


186 


80 


Ridge 2 


22:35:59.85 


+33:58:04.00 


318 


249 




Ridge 3 


22:35:59.73 


+33:57:37.16 


86 


41 




Bridge'' 


22:36:01.56 


+33:58:23.30 


175 


57 





Note. — Pointing positions and useful integration times of the CO spectra 
presented in this paper. The CO(l-O) and (2-1) lines were observed with the 
TRAM 30m telescope, and the CO(3-2) with the APEX telescope. 
" Starburst region at the northern tip of the SQ ridge 
Bridge feature that connects the ridge to the AGN NGC 7319 



3. CO OBSERVATIONS 
3.1. CO(l-O) and (2-1) with the IRAM 30m telescope 

The observations were made on June 28-30, 2009 at the 
IRAM 30 m telescope at Pico Veleta, Spain. We used the 
heterodyne receiver, EMIR commissioned in March-April 
2009. EMIR provides a bandwidth of 4 GHz in each of the two 
orthogonal linear polarizations for the 3,2, 1.3 and 0.9 mm at- 
mospheric windows. The four EMIR bands are designated as 
E090, E150, E230, and E330 according to their approximate 
center frequencies in GHz. Given the wide range of CO veloc- 
ities (6000-6700 kms"') exhibited by previous observations, 
and the broad linewidth of the S(l) H2 line in the shock, this 
increase in bandwith was crucial to our study. 

We observed the CO(l-O) and CO(2-l) Hues at 11 5. 27 GHz 
and 230.54 GHz by connecting the third and fourth parts of 
the WILMA backend to the horizontal and vertical polariza- 
tions of the E230 band. We focus this paper on the CO(l-O) 
results, since we did not sample the (1-0) beam with the (2-1) 
observations. To increase the redundancy of the CO(l-O) data, 
we used in parallel the 4 MHz filterbank with both polariza- 
tions connected to the E090 band. The receivers were tuned 
to 112.86 GHz and 225.72 GHz, which corresponds to a re- 
cession velocity of 6400 km s~' . The observations were done 
in wobbler switching mode, with a wobbler throw of 120" in 
azimuthal direction. This observation mode provides accurate 
background subtraction and good quality baselines, which is 
critical for measuring broad lines. Pointing and focus were 
monitored every 2h in stable conditions and every Ih during 
sunrise. The pointing accuracy was ~ 3" and the system tem- 
peratures were « 150-200 K at 115 GHz on the scale. 
The forward efficiency of the telescope was 0.94 and 0.91 at 
115 and 230 GHz and the main beam efficiency was 0.78 and 
0.58, respectively. The half-power beam size is 22" at 115 
GHz and 11" at 230 GHz. 

The data was reduced with the IRAM GILDAS- 
CLASS90'^^ software package (version of December 2010). 
The lowest quality spectra (typically with an RMS noise 
cr > 30 mK for a 3 min on-source integration) were rejected. 
We subtracted linear baselines. The final spectra presented in 
this paper are smoothed (with a banning window) to a veloc- 
ity resolution of 40 km s"' . We experienced problems with 
the vertical polarization. These spectra, which showed spuri- 
ous features, ripples and curved baselines, have been rejected 

For details about the Eight Mixer Receiver instrument, see 
http://www.iram.es/IRAMES/mainWiki/EmirforAstronomers 

Grenoble Image and Line Data Analysis System, 
http://www.iram .f r/l RAIVl F R/G I LDAS/ 



for this analysis. 

We observed five different positions in the center of the SQ 
group. The blue circles on Figure 1 show the CO(l-O) half- 
power beams (22" in diameter) overlaid on a HST WFPC2 
image of SQ. The three positions in the ridge are designated 
as Rl, R2 and R3, from North to South. The two other posi- 
tions are centred on the SQ-A star forming region and on the 
middle of the H2-bridge feature in between NGC 7319 and 
the main shock. The red lines show the contours of the H2 
0-0S(3) Hne emission detected by Spitzer (Cluver et al. 2010). 
Note that the half power beam width at 1 15 GHz matches very 
well the width of the H2 emission in the ridge. Table 1 lists 
the coordinates of these positions and the total ON+OFF time 
spent at each position. A total of w 20 hours has been spent 
for ON+OFF observations towards SQ. 

3.2. CO(3-2) with the APEX telescope 

Observations towards the Stephan's Quintet in the CO(3- 
2) transition were carried out during night time in July 2010 
with the Atacama Pathfinder Experiment telescope (APEX, 
Glisten et al. 2006), using the First Light APEX Submillime- 
ter Heterodyne receiver (developed by Max Planck Institut far 
Radioastronomie, MPIfR, and commissioned in May 2010) at 
345 GHz, FLASH345''*, along with with two newly commis- 
sioned MPIfR X-Fast Fourier Transform Spectrometer back- 
ends, XFFTS (Klein et al., in preparation), of 2.5 GHz band- 
pass each, and whose combination covers a 4 GHz bandpass. 
The beam of the telescope is 18.1" at 345 GHz (Glisten et al. 
2006). Focus was checked at the beginning of each night, 
and pointing was checked on planets (Jupiter, Saturn, and/or 
Uranus), as well as on CRL2688 every 90 minutes at least. 

System temperatures were in the range 180-270 K and 180- 
230 K in the respective cases of the Ridge 1 and SQ-A posi- 
tions, on the T^ scale. Antenna temperatures were converted 
to main beam temperatures using forward and beam efficien- 
cies, respectively, of 0.95 and 0.67 at 345 GHz. We used 
32768 channels for each XFFTS backends, yielding an effec- 
tive spectral resolution of about 76.3 kHz. The observations 
were performed in wobbler switching mode, with a throw of 
120" in the azimuthal direction. 

Again, the data was reduced by means of the IRAM 
GILDAS software package. We noted that w 15 channels (out 
of a total of 32768 covering the 4800-8400 km s"' velocity 
range) were affected by spurious spikes, likely due to the elec- 
tronics of other receivers that were powered-on at the time of 
the observations). Therefore, a sigma-clipping has been ap- 
plied. By scanning across all spectra, the sigma associated 
with each channel, (7^, is determined. The data is flagged if 
the signal is greater than 5ac and greater than 3.5 (t,, where 
(Tj is the sigma associated with one spectrum, since the signal 
does not reach this level. 

The baseline subtraction was always linear, as we enjoyed 
very good baseline stability due to excellent weather condi- 
tions during the observations. The final spectra shown here 
are smoothed to a velocity resolution of 40 km s"', for com- 
parisons purposes with the other transitions obtained at the 
IRAM 30m-telescope. 

4. KINEMATICS AND EXCITATION OF THE CO GAS 
4.1. CO spectra and analysis 

In the 345 GHz band the dual-polarization receiver FLASH operates a 
2SB SIS mixer provided by IRAM (Maier et al. 2005). 
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Table 2 

CO(l-O) observations with tlie IRAM 30m telescope: results from Gaussian decompositions 



Target region 


vco(i-O) 
[km s"^] 


^■vco(i-O) 
[km s^i] 


2^CO(l-0) 
\K km li^h 


M(H2) " 

[XiU ivi0j 


Lelgsj 


M(H^)/M(H;;)" 


SQ-A 


5987 ±14 
6681 ±22 
6891 ±11 


69±8 
230 ±35 
104 ±12 


0.9 ±0.1 
1.7±0.3 
1.6±0.3 


3.1±0.5 
5.7± 1.1 
5.3±0.9 


7.8±3.1 £54 
1.6±0.8£56 
3.1 ±1.2 £55 




SQ-A (sum) " 






4.2 ±0.4 


14.1 ±1.5 


2.0±0.8£56 


0.9 ±0.2 


Ridge 1 


t^niA _i_ 1 1 
6034 ±11 

6288 ±22 

6581 ±20 

6727 ±15 


144 ± 10 
700 ± 60 
71±8 
118±13 


1.9 ± 0.2 
1.5±0.2 
0.15±0.02 
0.4 ±0.06 


6.3 ± 0.8 
5.0 ±0.7 
0.50 ±0.08 
1.2±0.2 


"7 1 _1_ 1 Q 17 
/.i it i.O £ JJ 

1.3 ±0.4 £57 
1.3±0.5£54 

9.4 ±0.4 £54 




Ridge 1 (tot) " 






4.0 ±0.3 


13.0± 1.4 


1.4 ±0.4 £57 


0.26 ±0.03 


Ridge 2 


5970 ±28 
6038 ± 12 
6637 ± 35 


377 ±36 
93± 10 
342 ± 28 


2.3±0.3 
0.56 ±0.08 
1.0±0.1 


7.5± 1.1 
1.9 ±0.3 
3.4 ±0.4 


5.8±1.9£56 
8.6 ±3.2 £54 
2.1 ±0.6 £56 




Ridge 2 (tot) " 






3.87±0.32 


12.8± 1.2 


8.0±2.0£56 


0.33 ±0.07 


Ridge 3 


5777 ±25 


276 ±19 


3.5 ±0.4 


11.6±1.4 


4.8 ±1.2 £56 


0.10 ±0.02 


Bridge 


6239 ±21 
6605 ± 20 


423 ± 24 
298 ± 19 


2.5±0.3 
2.6 ±0.4 


8.2 ±0.9 
8.5 ± 1.0 


7.9±1.8£56 
4.1 ±1.0 £56 




Bridge (tot) " 






5.1±0.5 


16.7 ± 1.4 


1.2±0.2£57 


0.55 ±0.09 



Note. — Observational results of the IRAM 30m EMIR CO(l-O) observations. This table gathers the results from the fitting 
of the CO(l-O) spectra in multiple Gaussian velocity components (see Figure 2). The first three columns list the central velocity, 
FWHM, and intensity of the Gaussian components. Results are given for each velocity components, as well as for the sum of the 
components. The velocity component at 6900 km s"' in SQ-A is a newly discovered feature (see text for details). The aperture 
of the beam is 1 . 1 3 X (22 arcsec) ^ = 547 arcsec - . 



" Ht gas masses calculated using the Galactic conversion factor of A'(H 2)//co = 2 x 10^" cm^^ [K km s"']"' 
Kinetic energy in random motions, calculated as £ = 3/2M(H2)it^q (see Sect. 6.1). 

Ratio between the H 2 mass derived from the Spitzer spectral mapping of the mid-IR H 2 lines in the group to the H 2 mass 
derived from the CO(l-O) hne intensity at the different positions in the IGM of SQ. 
Sum over all velocity components 



Table 3 

CO(l-O) observations with the IRAM 30m telescope: results from integrations of the spectra over velocity ranges 



Target region 


Velocity range 

Vmin-Vmax [km S"'] 


(v) 
[kms-'] 


1cO{\-0) 

[Kkms-'] 


M(H2) " 
[xlO** Mq] 


[ergs] 


M(Hp/M(H5)" 


SQ-A 


5900-6200 
6200-6500 
6500-6800 
6800-7100 


5984 
6375 
6676 
6882 


0.9±0.1 
0.4 ±0.05 
1.5±0.2 
1.7±0.2 


3.0 ±0.2 
1.3±0.9 
4.8 ±0.5 
5.6±0.6 


3.2±0.8£54 
3.5±0.8£55 
7.4± 1.6£55 
1.7±0.5 £55 




SQ-A (sum) 






4.5 ±0.3 


14.7 ± 1.2 


1.3±0.2£56 


0.8 ±0.2 


Ridge 1 


5500-5900 
5900-6200 
6200-6500 
6500-6800 
6800-7500 


5831 
6041 
6325 
6656 
7039 


0.4 ±0.05 
2.5 ±0.3 
0.8±0.1 
0.9±0.1 
0.5 ±0.05 


1.3±0.3 
8.3± 1.0 
2.5±0.1 
2.9±0.1 
1.6±0.2 


2.2 ±0.6 £55 
1.1±0.4£56 
5.5±0.3 £55 
6.4 ±0.3 £55 
1.0 ±0.4 £56 




Ridge 1 (sum) 






5.2 ±0.4 


16.6± 1.5 


3.5±0.8£56 


0.20 ±0.03 


Ridge 2 


5700-6200 
6200-6500 
6500-6800 


5971 
6346 
6660 


2.6 ±0.3 
0.4 ±0.05 
0.7±0.1 


8.6± 1.2 
1.5±0.3 
2.4 ±0.3 


3.5±0.8£56 
3.8±0.9£55 
4.4 ± 1.0 £55 




Ridge 2 (sum) " 






3.7 ±0.4 


12.5± 1.5 


4.3± 1.1 £56 


0.34 ± 0.04 


Ridge 3 


5500-6100 


5788 


3.5 ±0.4 


11.6± 1.7 


4.5 ±0.6 £56 


0.10 ±0.02 


Bridge 


5900-7000 


6435 


4.9 ±0.4 


16.0 ±2.0 


2.4 ±0.5 £57 


0.57 ±0.09 



Note. — Idem as Table 2 but for integration of the CO(l-O) spectra over velocity ranges. The first three columns list the 
velocity range, the first moment velocity (mean, intensity-weighted velocity), and CO(l-O) intensity (cf. Eq. 1). Results are given 
for each velocity ranges, as well as for the sum of the integrations. 

" H2 gas masses calculated using the Galactic conversion factor of ^(H 2)//co = 2 X 10-" cm^^ [K km s"' 
Bulk kinetic energy, calculated from Eq. 3 (see Sect. 6.1). 

Ratio between the H 2 mass derived from the Spitzer spectral mapping of the mid-IR H 2 lines in the group to the H 2 mass 
derived from the CO(l-O) hne intensity at the different positions in the IGM of SQ. 
Sum over all velocity ranges at a given position. 
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Figure 1. CO and H2 observations in the Stephan's Quintet group. The grey-scale background image is a HST WFC2 V-band image. The blue circles (22 
ai'csec in diameter) show the IRAM 3()m CO(l-O) half-power beams where the observations presented in this paper have been taken. The corresponding spectra 
are shown on Figure 2. CO emission is detected at all these positions. The intruder galaxy, NGC 7318b, is colhding with intra-group gas (tidal debris from 
NGC 7319) at a relative velocity of ~ 900km s"'. The CO(l-O) green contours from BIMA observations by Gao & Xu (2000) correspond to levels of 2.4, 3.0, 
4.0, 5.0, 6.0, 7.0, and Sc (Icr = 0.7 Jy km/s). No CO was detected in the ridge because of limited sensitivity and bandwidth. Also note that the two ~ 3 -4(t CO 
features in the south were just outside (but near the edge) of the BIMA piimary beam. The red contours are the H 1 9.7 /im S(3) line contours from Spitzer IRS 
spectral mapping by Cluver et al. (2010), which clearly shows the kpc-scale shock produced by the galaxy collision. The H 2 contours are [1.0, 1.6, 2.3, 2.9, 3.6, 
4.2, 4.9]xl0"'' Wm-2 sr"'. 
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Figure 2. CO(l-O) and CO(2-l) spectra (left and right columns) observed with the IRAM 30m telescope at the five positions observed in the IGM of the SQ 
group (see Figure 1). CO emission with multiple velocity components is detected at all positions: in the shock, in SQ-A, and in the bridge regions. The red 
hne shows the result of the line fitting (see Tables 2 and 4 for the decomposition in multiple Gaussian components for each transition). The spectra have been 
smoothed to a spectral resolution of 40 km s"' . Note that the CO(l-O) and CO(2-l) beams do not match, and that the detection of the CO(2-l) line in the Ridge 3 
position is tentative (in this case the Gaussian fitting -dashed line- is not reliable). The aiTows indicate the pre-shock Hi velocities for the intruder galaxy and the 
intra-group medium. 
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Table 4 

CO(2-l) observations with the IRAM 30m telescope: results from Gaussian 
decompositions 



Target region ^C0(2 1) 
[km s J 


AvcO(2-l) 

rl -In 

[km s '] 


2cO(2-l) 
[Kkms"'] 


SQ-A 


5966 ± 25 
6719 ±36 
6906 ± 14 


81±8 
292 ± 35 
62 ±12 


0.76±0.11 

2.3 ± 0.4 

1.4 ±0.2 


SQ-A (tot) 






4.5±0.3 


Ridge 1 


5783 ±14 
6041 ± 16 
6101 ±39 
6721 ±12 


80±10 
92 ±12 
367 ± 34 
105 ±13 


0.50 ±0.08 
1.9±0.3 
2.4 ±0.4 

0.74±0.12 


Ridge 1 (tot) 




5.5±0.5 


Ridge 2 


5972 ±28 


252 ±36 


2.5 ±0.5 


Ridge 3 


5743 ± 25 


117±19 


1.9±0.4 


Biidge 


6318 ±22 
6552 ±46 
6606 ± 24 


77 ±12 
648 ± 36 
91±11 


0.7 ±0.1 
5.2± 1.6 
0.8±0.1 


Biidge (tot) 






6.7 ±1.6 


Note. — Idem as Table 2 but for the IRAM 30m EMIR CO(2-l) observa- 
tions. The size of the CO(2-l) beam is 1.13 X (11 arcsec) - = 137 arcsec^. 




Table 5 

CO(3-2) observations with APEX: results 




Target 


VCO(3-2) 
[kms"'] 


AvcO(3-2) 

[kms"'] 


2cO(3-2) 
[Kkms"'] 


SQ-A ■> 


6858 ±66 


80±31 


0.32±0.12 


Ridge 1 


5976 ± 15 


134 ±43 


0.65±0.13 


Note. 


Idem as Table 2 and 4 but for the APEX CO(3-2) observations 



(see Figure 3) obtained with the FLASH345/XFFTS receivers. The aperture 
of the CO(3-2) beam is 1.13 X (18.1 arcsec) ^ 370 arcsec -. 
" Tentative detection at a 2.3 cr level. 
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Figure 3. CO(3-2) spectra observed with APEX at the position Ridge 1 in 
the IGM of the SQ group (see Figure I). The red line shows the result of 
the line fitting (see Tables 5). The spectrum has been smoothed to a spectral 
resolution of 40 km s"' . 

In Figure 2 we show the CO(l-O) and CO(2-l) spectra for 
each of the five observed positions. The spectra are plotted 
over a 4000 km s"' bandwidth (out of the « 8000 km s"' 
bandwidth available at 115 GHz) at a spectral resolution of 
40 km s"' . The noise levels in the average spectra are typically 
« 0.7-1.6 mK for CO(l-O), and w 2-5 mK for CO(2-l) at 
this resolution. The two transitions of CO are detected with a 
high signal-to-noise ratio for all the targeted regions, except a 
2.5 (T tentative detection of CO(2-l) for the R3 position. Such 
detections are in sharp contrast with the previous interfero- 




5000 5500 6000 6500 7000 7500 8000 

Velocity [km s"'] 



Figure 4. Overlay of the CO(l-O), CO(2-l) and CO(3-2) spectra observed in 
the Stephan's Quintet shock, at the position Ridge 1. No beam size convo- 
lution has been apphed. The spectra have been shifted vertically for clarity. 
The arrows indicate the pre-shock Hi velocities for the intruder galaxy and 
the intra-group medium. 

metric observations, where no CO was detected in the ridge, 
probably because of limited velocity coverage and sensitivity. 
It shows that CO gas is present in the H2-luminous shock, and 
is co-spatial with the X-ray emitting hot plasma in this region, 
at a ~ 5kpc scale. 

The spectrum resulting from our pilot study with APEX at 
the ridge 1 position, where the CO(3-2) line is detected, is 
shown in Figure 3. In SQ-A, the line is tentatively detected at 
a 2.3 CT significance, and we consider this as an upper limit on 
the CO(3-2) flux. 

To compute the CO line parameters, we applied two differ- 
ent methods. First, we decomposed the CO spectra in multiple 
Gaussian components using the IDL non-linear least-squares 
MPFIT routine (Markwardt 2009). The red lines on Figures 2 
and 3 show the fitting result, which is used to derive the cen- 
tral velocities of the CO components, the line widths and in- 
tegrated intensities 2co- We checked that the residuals are 
below the 2 cr noise level for all positions and all lines. The re- 
sults of these Gaussian decompositions are given in Tables 2, 
4 and 5. 

Second, the CO line intensities Ico were derived from the 
integration of the CO main beam temperature, Tmb [K], over 
velocity ranges, from Vmi,, to v^ax, defined for each spectrum 
in Table 3, as the following: 

Ico[Kkms-']= ^ W^bdw, (1) 

J (v>-v„,i„ 

where (v) is the first moment (mean, intensity-weighted) ve- 
locity over that velocity range. This method allows us to better 
take into account the asymmetry in the CO profiles. The re- 
sults of these integrations are given in Table 3 and discussed 
in the following sections. 

4.2. Multiple CO velocity components 

The spectra in Figures 2 and 3 clearly show that multi- 
ple, broad velocity components are detected, pointing out the 
complexity of the kinematics of the CO gas in the SQ group. 
The CO line emission extends over a very wide velocity range 
w 5700-7000 km s"' (and perhaps until 7400 km s"' for the 
ridge 1 position). The four main velocity components are at 
5700, 6000, 6700 and 6900 km s"' (see table 2, 4 and 5 for the 
central line velocities derived from line fitting). The gas asso- 
ciated with the intruder NGC 7318b (5700-6000 km s"') and 
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with the intra-group medium (6700 km s"') are both detected 
in the ridge and in the eastern feature (the bridge) towards 
NGC 7319. These two CO velocity systems match with the 
velocities of the Hi gas observed in the group (Sulentic et al. 
2001; WilHams et al. 2002). In Ridge 3, the southern part of 
the shock, most of the gas seems to be associated with the in- 
truder, since the 6700 km s~' component is not detected. The 
velocity component at 6900 km s"' is a new feature that we 
detect only in the SQ-A region. It was not seen in previous CO 
or Hi observations because of limited velocity coverage, and 
its origin is still unknown. It could be the result of dynamical 
effects during the previous encounter between NGC 7319 and 
7320c, but current hydrodynamical simulations cannot repro- 
duce this feature. We note that an Ha velocity component at 
7000 km s"' was reported by Sulentic et al. (2001) in a clump 
to the North-East of SQ-A, but it is outside of our CO(l-O) 
beam. 

Gas at intermediate velocities, in between that of the in- 
truder and that of the intra-group, is detected. A weak, broad 
component centered at ss 6400 km s"' is detected in the 
ridge 1, ridge 2 and in the bridge positions. Interestingly, this 
intermediate component seems to be absent (or very weak, at 
a w 2ct significance) in SQ-A and vanishes at the southern 
end of the main shock region (ridge 3). The 6400 km s"' CO 
velocity component is consistent with the central velocity of 
the broad, resolved H2 S(l) line at 6360 ± 100 km s"' (see 
Appleton et al. 2006; Cluver et al. 2010, for an analysis of the 
high-resolution Spitzer IRS spectrum). This is the first time 
CO gas is detected at these intermediate velocities, and this 
gas is spatially associated with the warm H2 seen by Spitzer. 
Xu et al. (2003); Iglesias-Paramo et al. (2012) reported that 
the low velocity component (intruder) of the optical emission 
lines in the shock and in the bridge is centered at intermediate 
velocities (between 5900 and 6500 km s"'). In the northern 
(SQ-A) and southern star-forming regions, at both ends of the 
large-scale shock, the low velocity component of the optical 
lines is consistent with the intruder velocity. 

In Figure 4 we compare the three CO line profiles at the 
ridge 1 and SQ-A positions. At both positions, the central 
velocity of the CO(3-2) line coincide with the brightest CO(l- 
0) and (2-1) components. The FWHM of the CO(3-2) line 
matches those of the brighter CO( 1-0) and (2-1) lines. 

4.3. High CO velocity dispersions 

In the ridge, the FWHM of the main CO components are 
of the order of 100-400 km s"' (see tables 2, 4 and 5). 
In the ridge 1 position (and perhaps ridge 2), the under- 
lying intermediate component at 6400 km s"' is extremely 
broad, « 1000 km s~', with possibly a hint of emission be- 
tween 6900 and 7400 km s"'. This CO velocity coverage of 
« 1000 km s"' is comparable to the relative velocity between 
the intruder and the IGM, and to the broadening of the H2 
S(l) Une (FWHM« 870 km s"') observed by Appleton et al. 
(2006); Cluver et al. (2010). 

In SQ-A, the linewidths are of the order of 100 km s"', with 
the smallest dispersion for the 6000 km s"' line. The broad- 
est component at 6700 km s"' seems to have a left wing ex- 
tending towards 6400 km s"', suggesting again that some of 
the CO gas belonging to the IGM has been decelerated in the 
shock induced by the intruder galaxy. Interestingly, for the 
three transitions, the CO linewidths are significantly smaller 
in SQ-A than in the ridge and bridge, and this is perhaps a 



clue to understand why there is much more star formation in 
this region than in the ridge (see discussion in sect. 7. 1). This 
observational picture fits within the interpretation of two col- 
liding gas flows, where the shear velocities between the flows 
are maximum in the central region of the contact discontinuity 
(ridge), and smaller at the edges of the main shock structure 
(see Figure 5 of Lee et al. (1996) and Guiflard, R (2010)). 

The spectrum in the bridge is perhaps the most striking. 
It shows a very broad signal in between the extreme veloci- 
ties detected in SQ-A (6000 and 6900 km s"'), see Figure 2. 
The CO(l-O) profile is well fitted by a sum of 2 Gaussians, 
whereas in the CO(2-l) spectrum, narrower velocity compo- 
nents (centered at 6300 and 6600 km s"' ) emerge on top of a 
broader (FWHMw 650 km s"') emission. 

4.4. CO excitation 

Since the beams of the CO(l-O), (2-1) and (3-2) observa- 
tions do not match, and given that the angular size of the CO 
emitting regions are not constrained, we cannot infer useful 
information on the CO excitation (from the CO line ratios). 
We generafly find a CO(2-l) to CO(l-O) line ratio of w 1, 
which suggests that the CO is not far from being optically 
thick. This has to be confirmed with a complete CO(2-l) map- 
ping of SQ. 

5. DISTRIBUTION AND MASS OF MOLECULAR GAS 
5.1. Warm H2 masses from the Spitzer mid-IR spectroscopy 

To derive the warm (T > 100 K) H2 masses from the Spitzer 
observations, we extracted IRS spectra within the areas ob- 
served in CO(l-O) with the CUBISM software (Smith et al. 
2007a). The H2 line fluxes, listed in Table 6, were derived 
from Gaussian fitting with the PAHFIT IDL tool (Smith et al. 
2007b). The H2 physical parameters are estimated with 
the method described in Guillard et al. (2009). The spec- 
tral energy distribution of the H2 line emission is modelled 
with magnetic shocks, using the MHD code described in 
Flower & Pineau des Forets (2010). The gas is heated to a 
range of post-shock temperatures that depend on the shock 
velocity, the pre-shock density, and the intensity of the trans- 
verse magnetic field. We use a grid of shock models (vary- 
ing shock speeds) similar to that described in Guillard et al. 
(2009), to constrain the physical conditions (density, shock 
velocity) needed to fit the H2 line fluxes. These fits are not 
unique, and in principle, the observed H2 excitation arises 
from a continuous distribution of densities and shock veloci- 
ties. We find a best fit for a pre-shock density hh = 10^^ cm"^*. 
The initial ortho-to-para ratio is set to 3, and the intensity 
of the pre-shock magnetic field is set to Bq = y/nn w 30 /iG, 
in agreement with Zeeman effect observations of Galactic 
molecular clouds (Crutcher 1999), and the value inferred from 
radio observations of the synchrotron emission in the SQ ridge 
(Xu et al. 2003; O' Sullivan et al. 2009a). The H2 line fluxes 
and the warm H2 masses are computed when the post-shock 
gas has cooled down to a temperature of 100 K. This temper- 
ature is chosen because below « 100 K, the gas does not con- 
tribute significantly to the rotational H2 emission. The warm 
H2 masses derived from the shock models depend on this tem- 
perature, the lower the temperature, the longer the cooling 
time, hence the larger H2 mass. However, the relative dif- 
ference in cooling times (hence in H2 masses) between 100 K 
and 150 K is < 20%, so this choice of temperature does not 
impact much the results given in Table 7. 

A combination of two shocks with distinct velocities is re- 
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Table 6 

H2 rotational line intensities 



Target region Ht S(0) 


H2 S(l) 


H2 


S(2) 


H2 


S(3) 


H2 


S(4) 


H2 S(5) 


SQ-A 


1.47±0.11 


6.26±0.13 
















ridge 1 


0.79 ±0.04 


11.0± 1.56 


3.19 


±0.25 


11.5 


±1.31 


1.97 


±0.27 


4.97 ±0.41 


Ridge 2 


0.68 ±0.04 


7.06 ±0.42 


2.75 


±0.11 


13.3 


±2.35 


1.78 


±0.33 


4.31±0.31 


Ridge 3 


0.49 ±0.07 


9.55± 1.51 


3.42 


±0.31 


5.52 


±1.53 


0.94 


±0.34 


2.85 ±0.39 


Bridge 


0.88 ±0.07 


4.88 ±0.08 

















Note. — H2 line fluxes derived from the Spitzer IRS spectra extracted at the positions of our CO 
observations. Fluxes are in units of 10"" W m"^. 



Table 7 

H2 line flux fitting with MHD shock models: results 



Target region 


shock velocities 
[kms-'] 


Mass flow 
[Mq yr-1] 


Cooling time 

[yr] 


M(H™™) 

[Mq] 




SQ-A 


5 


1.47 X lO'' 


8.09 X lO-* 


1.19 X 10' 


Ridge 1 


8 

35 


6.28 X 10^ 
5.45 X 10^ 


5.33 X lO* 
9.77 X 10^ 


3.35 X 10** 
5.33 X 10* 


Ridge 2 


6 


6.61 X 10^ 


6.30 X lO* 


4.16 X 10** 


35 


5.40 X 10^ 


9.77 X 10^ 


5.27 X 10* 


Ridge 3 


13 


3.32 X 10^ 


5.52 X lO-* 


1.17 X 10** 


40 


1.74 X 10- 


7.76 X 10^ 


1.35 X 10* 


Biidge 


5 


1.13 X lO"* 


8.09 X 10* 


9.13 X 10** 



Note. — H2 fine fluxes are fitted with low-velocity MHD shocks in dense 
molecular gas, for a pre-shock density of rm = lO'cm"' and a pre-shock 
inagnetic field of 30 ^G. The table gives the best-fit shock velocities, mass 
flows, cooling times and warm H2 masses. 

quired to match the observed H2 line fluxes. As discussed in 
Guillard et al. (2009), the H2 masses are derived by multiply- 
ing the gas cooling time (down to 100 K) by the gas mass flow 
(the mass of gas swept by the shock per unit time) required to 
match the H2 line fluxes. The shock model parameters, gas 
cooling times, mass flows, and warm H2 masses are quoted in 
Table 7. The IRS spectra in different regions of the shock, and 
a detailed description of the spatial variations of H2 excitation 
(full 2-dimension excitation diagrams) in the SQ intra-group 
medium will be given in a separate paper (Apple ton et al., in 
preparation). 

5.2. H2 masses derived from CO observations 

We derive the masses of H2 gas from the CO(l-O) 
line observations, using the Galactic conversion factor X = 
Af(H2)/Xco = 2 X 10-" cm-2 (K km s"')-': 

A?(H2) 

^co " 



Ml 



H2 



Ico 



xd xftxlmn 



(2) 



7-obs 



Mil-, » 
—=3.3x10**- 
M0 Kkms" 



A^(H2)/Ico / d 



2x102" V94Mpc 



Php 
22" 



where I^q is the observed velocity integrated CO(l-O) line 
intensity, d is the distance to the source, and il is the area cov- 
ered by the observations in arcsec", with = 1.13 x for 
a single pointing with a Gaussian beam of Half Power beam 
width of To derive the total mass of molecular gas, one 
has to multiply Me, by a factor 1.36 to take the Helium con- 
tribution into account. Tables 2 and 3 give the H2 gas masses 
for the different positions, derived from the CO line intensi- 
ties. In the following, we discuss the masses of molecular 
gas detected in the three main areas, the ridge, SQ-A, and the 
bridge. 



In the SQ ridge (positions Rl, R2, and R3), summing 
over all the CO(l-O) velocity components, and assuming a 
Galactic CO to H2 conversion factor, we derive a mass of 
« 4.1 X 10^ Mq of H2 gas, corresponding to a surface den- 
sity of i;(H2) = 12.0 ± 1.3 Mq pc"2. Note that these three 
pointings cover w 75% of the whole H2-bright ridge (the Rl 
and R2 beams are partially overlapping). 

Within the CO(l-O) aperture n= 1.13 x (22 arcsec)^ = 
547 arcsec^ centered on SQ-A, we find a total H^ mass of 
1.3 X 10*^ Mq. Lisenfeld et al. (2002) find w 3 x lO** Mq over 
a much larger area (60 x 80 arcsec^) that partially overlaps the 
north of the ridge. However, because of a limited bandwith, 
they only detected the 6000 and 6700 km s"' velocity compo- 
nents, and not the higher velocity component at 6900 km s"' 
that contains almost half of the mass in the core region of SQ- 
A (see Table 2). Based on the MIPS 24 /im map, we estimate 
the area of SQ-A to be sa 40 x 40 arcsec^. If we scale the mass 
we find in our aperture fl to this larger aperture, we find a total 
mass of H2 of 3.8 x 10^ Mq. This is comparable to the mass 
derived by Lisenfeld et al. (2002), because the larger aperture 
used to estimate their masses compensates their non-detection 
of the 6900 km s"' component. 

In NGC 7319's bridge, summing over all velocities, we find 
Mh, w 1.7 X 10^ Mq within the aperture O. Note that this 
aperture matches quite well the width of the bridge, but does 
not cover its whole extension towards the AGN. Therefore, 
the total mass of molecular gas in the bridge may be a fac- 
tor of « 2 larger Smith & Struck (2001) find « 5 x 10'' Mq 
within a half -power beam of 55 arcsec in diameter centered 
on the nucleus of NGC 7319. This beam overlaps our bridge 
region. The global morphology of the CO gas in SQ is still 
unknown, but it is surely linked to the geometry of the col- 
lision and the morphology of the colliding gas flows. For 
instance, the origin of the molecular gas in the bridge struc- 
ture is still an open question. It is most likely a result of a 
previous tidal interaction (with NGC 7320c for instance) that 
would have stripped some material from NGC 73 19's galactic 
disk. Then the interaction between the new intruder and this 
"tidal bridge" may trigger further molecular gas formation in 
this region, since the proportion of post-shock gas (with re- 
spect to the pre-shock) in this region is the highest observed. 
Aoki et al. (1996) also reported from optical spectroscopy that 
an outflow from the AGN is present, which seems oriented 
in the direction of the bridge feature. Some of the molecu- 
lar material present in the bridge may be gas entrained from 
NGC 7319's disk by the outflow. However, it is unlikely that 
the outflowing ionized gas can entrain the dense molecular 
gas on such large distances (> 10 kpc). 

These observations represent a substantial revision of the 
mass and energy budgets of the collision. Previous interfer- 
ometric observations may have missed the molecular gas in 
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the ridge because of limited sensitivity and bandpass and be- 
cause the CO emission in the ridge may be more diffuse than 
in SQ-A. We detect at least 4 times more CO emission in the 
SQ ridge and in the bridge than in the starburst region SQ- 
A. Although these observations do not cover all the IGM, the 
CO emission seems distributed over a 40 kpc scale along the 
ridge, and in total (SQ-A + ridge + bridge), > 5 x 10^ Mq of 
H2 gas is lying outside the optical disks of the group galax- 
ies. Note that these masses depend on the X-factor, which is 
poorly constrained (sect 5.3). 

5.3. Uncertainties in the H2 masses 

We computed the uncertainties in the warm H2 masses de- 
rived from the Spitzer spectroscopy by fitting shock models 
to the lower and upper values of the H2 line fluxes, and we 
found uncertainties on the H2 masses of the order of 20% 
(based solely on the observational uncertainties). The warm 
H2 masses derived from shock models are 10-35% larger 
than local thermal equilibrium LTE models (we performed 
single temperature fits for SQ-A and the Bridge regions, and 
two temperatures fits for the positions in the ridge). This is 
mainly because (i) the values of the H2 ortho-to-para ratio are 
different in the two models, (ii) at «h =10^ cm"^^, the S(0) and 
S(l) lines are not fully thermalized, and (iii) the H2 mass de- 
rived from shock models is dependent on the final temperature 
at which the H2 level populations are calculated. 

We derived H2 masses from the CO(l-O) line intensities us- 
ing the Galactic value of the CO-to-H2 mass conversion factor 
(the X-factor), because the metallicity in the IGM is slightly 
above the solar value (Xu et al. 2003). However, the X-factor 
is uncertain, essentially because it is poorly constrained in 
shocked environments (Downes & Solomon 1998; Gao et al. 
2003), and more generally in low density molecular gas. For 
instance, in a similar situation in the bridge between the Taffy 
Galaxies (UGC 12914 and 12915), Braine et al. (2003) esti- 
mated a X -factor roughly five times below the Galactic value 
using '^CO measurements. The dust emission modelling by 
Natale et al. (2010), using their SpitzerMWS 70 and 160 ^m 
measurements in the shock, gives S(H2) = 9.8 M0 pc~^ (as- 
suming a Galactic dust-to-gas mass ratio Za = 0.0075). This 
value is close (« 20 % smaller) to the average H2 mass sur- 
face density derived from the CO measurements in the ridge 
(I](H2) = 12.0 ± 1.3 M0 pc"^), assuming a Galactic conver- 
sion factor (Table 2). This is within the uncertainties of the de- 
rived values from the dust modeling. The large uncertainties 
on the far-lR photometry in the shock are dominated by the 
difficulty to spatially separate the far-lR emission of the ridge 
itself from the bright sources surrounding the shock. Upcom- 
ing Herschel observations will help in estimating the X-factor 
in the shock and in the surrounding star forming regions. We 
defer this discussion to a future paper The uncertainty on 
the X-factor affects the mass and energy budget discussed in 
Sect. 7.3 but does not change the general conclusions of this 
paper. 

6. ENERGETICS AND EXCITATION OF THE MOLECULAR GAS 

6.1. Kinetic energy carried by the molecular gas 

A first estimate of the kinetic energy of the molecular 
gas was based on Spitzer/IRS spectroscopy of the H2 line 
emission from SQ (Guillard et al. 2009). Assuming that the 
8.8 X 10^ M© of warm H2 gas detected by Spitzer has a ve- 
locity dispersion of 370 km s"' (Cluver et al. 2010), the tur- 
bulent kinetic energy carried by the warm H2 gas is w 3.6 x 



10^^ ergs. Given the limited spectral resolution of Spitzer/IRS 
(w 600 km s~'), this estimate was crude. The present CO 
observations allows us to compute the kinetic energy more 
accurately. 

As for the CO fluxes and masses, we estimate the kinetic en- 
ergy carried in the CO emitting gas with two methods. First, 
from the Gaussian decompositions, the turbulent kinetic en- 
ergy of the H2 gas is derived from = 3/2MH,cr^O' where 
Mh, is the H2 mass derived from the CO(l-O) flux (Eq. 2) and 
(Jco the velocity dispersion along the line of sight (measured 
on the Gaussian fits of the CO(l-O) line). The factor of 3 in 
^kin takes into account the three dimensions. These kinetic 
energies are listed in Table 2 for each of the CO velocity com- 
ponent and each observed positions. 

The second method consists in integrating Tmh x (v- (v))^ 
over velocity ranges, from Vmin to Vmax, defined for each spec- 
trum and each velocity components in Table 3: 



£kin=Tx3.3xlO'^M0 



-(v) 



N(H2)/X co 
2x10 



20 



94Mpc 



22" 



X / (v-(v))2x7;,bdv, 

'(v>-v„™ 



(3) 



where (v) is the mean, intensity-weighted (first moment) ve- 
locity. The comparison of Tables 2 and 3 shows that the two 
methods give similar results, essentially because we defined 
velocity ranges that bracket the main peaks of CO emission. 
Integrations over these velocity ranges give lower limits of the 
bulk kinetic energy of the H2 gas for each positions. An inte- 
gration over the full extent of the CO velocity range gives an 
upper limit (see Table 3). 

The CO emitting gas carries a huge amount of kinetic en- 
ergy, « 4 X 10^^ ergs when summing on all positions. The 
kinetic energy of the H2 gas is at least 5 times higher than the 
thermal energy of the hot plasma heated by the galaxy col- 
Usion (Guillard et al. 2009, 2010, estimated that this thermal 
energy is w 9 x 10^^ ergs for the entire ridge and bridge). Be- 
fore the collision, much of the kinetic energy of the gas is 
associated with the motions of the intruder and intra-group Hi 
gas. The pre-shock Hi mass in the area covered by the CO ob- 
servations is estimated to be 1.0-2.8 x 10^ Mq (based on an 
extrapolation of the Hi observations in NGC 7319's tidal tail 
to the area of the shock, see sect. 2.2. 1 of Guillard et al. 2009). 
Assuming that this gas is shocked at a velocity of 600 km s"' , 
the pre-shock Hi kinetic energy is0.5-1.3x 10^^ ergs. There- 
fore, the CO emitting gas is carrying a large fraction of the 
pre-shock kinetic energy of the gas. Most of this kinetic en- 
ergy has not been thermalized yet. We suggest that this turbu- 
lent kinetic energy is the energy reservoir that powers the H2 
line emission and maintains a large fraction of the mass of the 
H2 gas warm in SQ. 

6.2. Turbulent heating of the molecular gas 

To assess the physical state of the molecular gas, we first 
compute the ratio between the warm H2 mass derived from the 
Spitzer/IRS spectroscopy and the H2 mass derived from the 
CO(l-O) observations, which we call "mid-IR/CO H2 mass 
ratio" . These ratios are listed for each observed position in 
Table 2. On average, this ratio is much higher in SQ than in 
star-forming galaxies, where the mass of cold molecular gas 
(< 100 K, usually traced by CO) is 25-250 times greater 
than the warm (> 100 K) H2 mass detected in the mid-lR 
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Table 8 

Intensities of the PAHs complexes 



Target region 


7.7 


11. 3 /im 


12.6/im 


17 


SQ-A 








4.67 ±0.21 


Ridge 1 


1.54 ±1.20 


8.53 ±0.48 


2.23 ±0.22 


I.84±0.36 


Ridge 2 


7.60 ±1.09 


7.59 ±0.26 


3.01 ±0.26 


2.79 ±0.32 


Ridge 3 


17.2 ±1.44 


9.31 ±0.27 


1.80±0.37 


1.50 ±0.27 


Bridge 








0.96 ±0.26 



Note. — PAH intensities derived from Spitzer IRS spectra at the positions 
of our CO observations, and integrated over the CO(l-O) beam. Fluxes are in 
units of 10''^ W mr^ 



(Roussel et al. 2007). The H2 mass in the ridge derived from 
the CO(l-O) line emission is only a factor of w 5 larger than 
the warm H2 mass seen in the mid-lR domain (8.8 x 10*^ M©, 
see Table 7). Surprisingly, this H2 mass ratio is extreme in 
SQ-A, whereas in the southern star-forming region (ridge 3) it 
is the lowest, with a value (0.1) compatible with star forming 
galaxies. 

The mid-lR/CO H2 mass ratios in SQ are comparable to 
those observed in molecular clouds near the Galactic center 
(Rodriguez-Fernandez et al. 2001), or in the bridge between 
the Taffy galaxies (UGC 12914/5), where 0.1 < M[j"V^h' ^ 
0.7 (Peterson et al., submitted), or in H2-luminous radio 
galaxies (Ogle et al. 2010). In all these environments, it has 
been concluded that the molecular gas is heated by shocks 
rather than UV photons. In such environments, the CO emis- 
sion may trace not only cold molecular gas (« 10 - 100 K) but 
also warmer gas (« 100 to a few 10^* K, like the gas traced by 
mid-lR rotational H2 lines). 

The H2 luminosity-to-mass ratio, £(H2)/M(H2) measures 
the cooling rate (per unit mass) of the molecular gas 
(Guillard et al. 2009; Nesvadba et al. 2010). We compare this 
ratio to the turbulent heating rate per unit mass, written as 
Ft = j(Jj/tdis (Mac Low 1999). or is the turbulent velocity 
dispersion of the H2 gas, and fdis is the dissipation timescale, 
often written as fjis = I far, where / is the size of the region 
over which the velocity dispersion is measured. If the H2 
emission is powered by the dissipation of turbulent energy, 
then the total energy dissipation rate must be at least equal to 
the H2 luminosity (Guillard et al. 2009): 



3 4 CCHoJ 



2 / 



(4) 



M(H2) 

On average in the SQ IGM, we find /:(H2)/M(H2) = 0.10 ± 
0.02 Lq Mq'. At the spatial scales of a GMC (/ « lOOpc), 
this condition implies aj > 35 km s"'. This is an order of 
magnitude higher than the CO velocity dispersion observed 
in Galactic GMCs at these spatial scales (Heyer et al. 2009), 
which shows that the molecular gas in SQ is much more tur- 
bulent than in the Galaxy. 

7. DISCUSSION 

These CO observations raise questions related to the effi- 
ciency of star formation in the IGM and the origin of the CO 
gas. 

7.1. Shocks, star formation, and physical state of the 
molecular gas in the SQ group 

Following Nesvadba et al. (2010), Figure 5 presents a di- 
agram similar to the classical Schmidt-Kennicutt relation- 
ship. Based on the empirical correlation between PAH emis- 
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Figure 5. Star-formation rate (assumed to be traced by the surface luminosity 
of the PAH emission) as a function of the surface mass of Hi, derived from 
CO(l-O) line measurements, integrated over the full velocity range. The solid 
hne shows the relationship obtained for the SINGS sample. Empty diamonds 
and filled dots mark dwarf galaxies, star-forming galaxies, and AGN. The 
SINGS data is from Roussel et al. (2007). The filled red pentagons show 
the positions within the Stephan's Quintet group where we performed CO 
observations and extracted IRS spectra. Rl, R2 and R3 are the three positions 
in the ridge. SQ-A denotes the northern star-forming region and Br the bridge 
(see Figure 1). For comparison, the orange triangles mark the H2-luminous 
radio galaxies of the Ogle et al. (2010) sample for which we have both clear 
PAH and CO detections. 



sion and star formation in starbursts (Calzetti et al. 2007; 
Pope et al. 2008), the 7.7 i^im PAH emission is used as a tracer 
of star formation. The intensities of the PAH complexes mea- 
sured on the IRS spectra are listed in Table 8. 

We find that the SQ shock and bridge regions are signifi- 
cantly offset (up to -1.9 dex) from the correlation obtained 
with the SlNGs sources (data taken from Roussel et al. 2007). 
Similarly low PAH-to-CO surface luminosity ratios are also 
found in some H2-luminous radio galaxies (Nesvadba et al. 
2010). We interpret this relative suppression of PAH emis- 
sion as evidence for a lower star formation efficiency. 

Part of this offset could come from an under-abundance of 
PAHs with respect to the Galactic value (possibly caused by 
their destruction in shocks). Our modelling of the dust emis- 
sion showed that the observed PAH and mid-lR continuum 
emission are consistent with the expected emission from the 
warm H2 gas in the shock for a Galactic dust-to-gas mass ratio 
and Galactic PAHs and very small grains (VSGs) abundances 
(Guillard et al. 2010). Since the H2 mass derived from the CO 
observations is only a factor of 5 (or less, given the uncertain- 
ties on the X-factor) higher than the H2 mass derived from IRS 
spectroscopy, it is unlikely that the PAH abundance would be 
depleted by more than an order of magnitude. Therefore, we 
favour the possibility that the efficiency of conversion from 
molecular gas to stars in the shock is much lower than in "nor- 
mal" star forming galaxies. On the contrary, the SQ-A region 
is consistent with these galaxies and with star-forming colli- 
sion debris in colliding galaxies (Boquien et al. 2010). In this 
region, the star formation rate is more than one order of mag- 
nitude higher than in the shock (Cluver et al. 2010). 

Interestingly, we find that the regions where the star for- 
mation is the most suppressed are those where the CO ve- 
locity dispersion is the highest. The very high velocity 
dispersion observed in the CO gas in the SQ shock and 
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Table 9 

Mass, energy and luminosity budgets of the pre-shock and post-shock gas in the SQ ridge. 





PRE-SHOCK GAS 




POST- 


SHOCK GAS 




gas phases 


Hot Plasma" 


Hi'' 




Hot Plasma" 


HII"^ 


Hi'' 


H2'= 


T[K] 
PAb [Kcm-3] 
Nh [cm--] 


3.2 X 10--^ 
5.7 X 10* 


3 X lO^o 


10- 10^ 


1.17 X 10-^ 
6.9 X lO* 
1.9 X 10^ 


1.4 X 10" 


< 5.8 X 10" 


10- lO"* 
5 X 102" 


Masses [Mq] 


2.6 X 10** 


1.0-2.8 X lO'' 


< 2.5 X 10'' 


5.3 X 10** 


4.9 X 10' 


< 2 X 10** 


1.6-4.1 X 10'' 


Energy [erg] 


Thermal 
3.4 X lO'f- 


Kinetic 
0.5-1.3 X lO^** 


< 1.9x 10" 


Thermal 
9 X lO^f- 


6 X 10^5 




Kinetic 
4-13 X 10^7 


Flux [W m-2] 


X-rays 
6.6 X 10-" 




CO(l-O) 
< 1.9x K)-" 


X-rays 
31.4 X 10-1^ 


Ha 
7.8 X 10-" 


Ol 

5.5 X 10-" 


mid-IR H: lines " 
7.5 X 10-"' 


£* [erg S-' ] 


6.2 X 10*' 




< 2 X 10*' 


2.95 X 10*' 


6.5 X 10"" 


4.6 X 10"" 


8 X 10"' 



Note. — All the numbers are scaled to the main shock aperture defined in Cluver et al. (2010): Ams = 77 x 30 arcsec^ = 37.2 X 13.2 kpc^. This aperture does 
not include the bridge region. The pre-shock gas is mainly in the form of Hi gas, contained in NGC 7319's tidal filament, and molecular gas. After the shock, the 
mass is distributed between the hot X-ray emitting plasma and the H2 gas (see discussion in Guillard et al. 2009). 
" Chandra observations of the extended X-ray emission in the shock and the tidal tail (O'Sullivan et al. 2009b). 
'' Based on extrapolation of Hi observations in the tidal tail by Sulentic et al. (2001); Williams et al. (2002). 
Ha and Ol optical line observations by Xu et al. (2003). 

From Spitzer IRS observations. The H2 line flux is summed over the S(0) to S(5) lines (Cluver et al. 2010). 
Derived from the CO observations presented in this paper, assuming a Galactic conversion factor 
^ Luminosities are indicated assuming a distance to SQ of 94 Mpc. 

bridge could maintain most of this gas in a diffuse and non- 
gravitationally bound state, which would explain why this 
molecular gas is inefficient at forming stars. This trend is 
also observed in other interacting systems like the tidal dwarf 
galaxy VCC 2062 (Due et al. 2007) or Arp 94 (Lisenfeld et al. 
2008), where the star formation is only observed in regions 
where the CO lines have a small velocity dispersion. Suppres- 
sion of star formation in highly turbulent molecular disks has 
also been observed in radio-galaxies (Nesvadba et al. 201 1). 

7.2. Origin of the CO gas in the IGM ofSQ: pre- or 
post-shock gas ? 

The CO gas in the IGIVI of SQ has two possible origins; it 
is either associated with the pre-shock or the post-shock gas 
(see discussion in Guillard et al. 2009). 

We call pre-shock molecular gas the H2 clouds that have not 
seen the fast (« 600 km s"' ) intercloud shock produced by the 
galaxy collision (materialized by the X-ray emitting ridge), 
or the clouds for which the fast intercloud shock have swept 
across them, but these clouds have survived to the shock. 
In the latter case, these pre-shock molecular clouds must be 
dense enough not to be ionized by the radiation field of the fast 
intercloud shock (propagating in the tenuous, «h ~ 10"^ cm""*, 
and X-ray emitting gas). Such clouds have too much inertia to 
be decelerated by the tenuous gas, and retain their pre-shock 
velocity (in Appendix we show that the deceleration timescale 
is longer than the age of the shock). Therefore, we expect to 
detect this CO gas at velocities close to the two pre-shock ve- 
locities (5700-6000 and 6700 km s"'). 

We call post-shock gas the molecular gas that is formed 
from diffuse gas (e.g. Hi gas) that has been heated by the 
transmitted shocks into the clouds, decelerated in the frame of 
reference of the intercloud shock front, and that had time to 
cool and become molecular within this deceleration timescale. 
This is the scenario that Guillard et al. (2009) proposed to 
explain the formation of H2 in the SQ ridge. This post- 
shock H2 gas is dynamically coupled to the lower density 
phases. In the frame of reference of the observer, the post- 
shock CO gas, formed from shocked intruder gas, is deceler- 
ated to velocities redder than the pre-shock intruder velocity 



(at V > 5700-6000 km s"'), while the post-shock CO gas, 
formed from shocked IGIVI gas'^, is decelerated to velocities 
bluer than the pre-shock IGM velocity (V < 6700 km s"'). 
Therefore, some post-shock CO gas should appear at a veloc- 
ity intermediate to that of the two colliding gas flows. How- 
ever, depending on the coupling efficiency and the dynamics 
of the galaxies in the group, some post-shock gas could also 
be found close to the IGIVI (NGC 7319's tidal debris) or the 
intruder (NGC 7318b) velocities. 

The detection of CO emitting gas at both pre-shock and in- 
termediate velocities between that of the intruder and that of 
the IGIVI shows that the molecular gas in the IGM of SQ orig- 
inates from both pre-shock clouds, and from the formation 
of molecular gas in the cooling post-shock gas that has been 
accelerated in the frame of the shock. The observed CO kine- 
matics show that the pre-shock components have lower ve- 
locity dispersions than the post-shock components (sect. 4.3), 
suggesting that the shock is driving turbulence into the post- 
shock gas. 

7.3. Pre- and post-shock gas masses and energy budgets 

The mass and energy budget of the pre-shock and post- 
shock gas in the ridge is summarized in Table 9, which is 
an updated version of Table 1 in Guillard et al. (2009), now 
taking into account the CO data, and the H2 spectral map- 
ping of the shock (ridge) region by Cluver et al. (2010). Dis- 
tinguishing spatially between the pre- and post-shock gas is 
difficult with the present CO observations. We rather rely 
on the kinematics to separate and quantify the masses and 
kinetic energies carried by each component. The uncertain- 
ties in the relative masses associated with both components 
are high, essentially because the dynamical coupling of the 
H2 gas with lower density phases and the global dynamics 
of the system are poorly constrained. In Table 9, the CO 
components with central velocities in the ranges 5700-6100 
and 6600-6800 km s"' are associated with the pre-shock gas 
of the intruder and the intra-group medium, according to Hi 

'^ CO should form on a similar timescale as H2 in the post-shock gas (see 
Figure C.2 of Guillard et al. 2009). 
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(Sulentic et al. 2001) and optical data (Xu et al. 2003). Out- 
side these ranges, the CO emission is counted as post-shock 
gas. These velocity boundaries are highly uncertain and are 
likely to underestimate the proportion of post-shock gas, and 
over-estimate the pre-shock gas masses. We thus consider 
the post-shock gas masses as lower limits, and pre-shock gas 
masses as upper limits. In particular we choose to include the 
6000 km s"' component in the pre-shock budget as suggested 
in Moles et al. (1997); Sulentic et al. (2001); Xu et al. (2003), 
but this is questioned by other authors (Williams et al. 2002; 
Lisenfeld et al. 2002) who argue this is intruder gas that has 
been accelerated by the shock from 5700 to 6000 km s"' . 

Under this assumption, the post-shock cold H2 mass is 
« 1/3 of the mass detected at the pre-shock velocities. Al- 
ternatively, if we consider that the 6000 km s"' component 
belongs to the post-shock gas, this proportion is w 3 /2. In ei- 
ther case, both pre-shock origin and post-shock CO formation 
are important. 

The kinetic energy associated with pre- or post-shock gas 
has been calculated from Eq. 3, based on the same velocity 
separation as for the masses (Sect. 7.3). For the post-shock 
gas, we quote the range of kinetic energy, being computed 
velocity component by component, or integrated over the full 
velocity range (after subtraction of the CO signal at the pre- 
shock velocities). 

8. CONCLUSIONS AND FINAL REMARKS 

We report on CO(l-O), (2-1) and (3-2) line observations 
with the IRAM 30m and APEX telescopes in the galaxy-wide 
shock of Stephan's Quintet, aimed at studying the CO gas 
counterpart to the powerful warm H2 emission detected by 
Spitzer/IRS. The main results are the following: 

• Multiple CO velocity components are detected, at the 
pre-shock velocities (namely that of the intruder galaxy 
at 5700 km s"' and that of the intra-group medium 
at 6700 km s"'), but also at intermediate (6000- 
6400 km s"' ), and higher (6900 km s"' in SQ-A) veloci- 
ties. The CO gas originates from both pre-shock molec- 
ular gas, and post-shock gas, i.e. diffuse gas heated by 
the shock that had time to cool and become molecu- 
lar within the collision age. CO lines are intrinsically 
broader in the ridge and in the bridge than in the SQ-A 
and southern star-forming regions. A large fraction of 
the gas kinetic energy of the collision is carried by the 
turbulent velocity dispersion of the CO emitting gas, 
which means that most of this energy has not been ther- 
malized yet. 

• The CO(l-O) and (2-1) emission is distributed over the 
intra-group shock (the ridge) and in an eastern bridge 
feature that connects the shock to NGC 7319. We also 
detected the CO(3-2) line in the ridge and tentatively in 
SQ-A. Large amounts of CO emitting gas are detected, 
« 5.4 X 10^ M0 (after conversion from CO(l-O) emis- 
sion to H2 using the Galactic value) in the ridge and 
in the bridge. Note that this mass could be lower by 
a factor of a few because of the large uncertainties on 
the CO to H2 conversion factor This gas is co-spatial 
with warm H2 seen in the mid-lR domain, H II gas, and 
X-ray emitting hot plasma. The ratio between the H2 
mass derived from the Spitzer spectroscopy to the H2 
mass derived from the CO observations is w 0.2, which 
is 1 - 2 orders of magnitude higher than in star-forming 



galaxies, showing that the heating rate of the molecular 
gas in SQ is particularly high. This high heating rate is 
likely to be powered by the turbulent dissipation of the 
kinetic energy involved in the galaxy collision. This in- 
terpretation implies that the velocity dispersion on the 
scale of giant molecular clouds (GMCs) in SQ is an or- 
der of magnitude larger than the Galactic value. 

• We find that the ratio of the PAH-to-CO surface lumi- 
nosities in the shock and in the bridge regions are much 
lower (up to a factor 75) than the observed values in 
star-forming galaxies. In the SQ-A star-forming region, 
the PAH-to-CO ratio is compatible with the classical 
Schmidt-Kennicutt law. We suggest that the star for- 
mation efficiency in the shock of SQ is lower than in 
"normal" galaxies. This could be a consequence of the 
high turbulence within the H2 gas. 

Our census of the molecular gas in the IGM of SQ is still 
incomplete. High angular resolution interferometric radio ob- 
servations of SQ will help to determine the dumpiness of the 
CO emitting gas and perhaps spatially separate the different 
velocity components. Far-infrared photometry is also needed 
to estimate the total dust mass, and thus the cold gas content 
independently of the CO to Ht conversion factor 

A global model for the formation and kinematics of the 
molecular gas present in the IGM of SQ can only be addressed 
by means of complex numerical simulations. A first step 
has been achieved by Renaud et al. (2010), who presented 
a pilot A^-body study with a solution that matches the over- 
all geometry of the group. Hydrodynamical simulations by 
Hwang et al. (201 1) took a step further including the gas com- 
ponent and made first headways in our understanding of the 
gas kinematics in the system. However, presently, the exist- 
ing numerical simulations do not include the molecular gas 
and cannot reproduce the CO and H2 observations of SQ pre- 
sented here and in Appleton et al. (2006); Cluver et al. (2010), 
where most of the molecular gas is not associated with star 
formation. The prescriptions used in the numerical models 
to describe the conversion from interstellar gas to stars will 
have to be refined, especially in the treatment of the gas cool- 
ing, energy exchanges between the phases of the interstellar 
medium, and turbulent cascade of the kinetic energy of the 
collision. 
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APPENDIX 

Dynamical coupling of the CO gas to other gas phases 

In the following we show that it is unlikely that the CO gas observed at intermediate velocities originate from pre-shock 
molecular clouds, essentially because the timescale of dynamical coupling, i.e. the time it takes for the pre-shock cloud to 
decelerate and become comoving with the shocked intercloud medium is much longer than the age of the collision. 

Let us consider the case of an H2 cloud of density pc and radius Rc, initially in pressure equilibrium with the inter-cloud 
medium (of density pi), that is run over by a fast shock wave of velocity V, in the inter-cloud medium. The timescale for 
the cloud compression due to the rise of the external pressure, so-called the crushing time, is given by fcmsh = Rc/yc, where 

Vc ~ \ff'^i ~ ^'1 velocity of the transmitted shock into the cloud (x = Pcj p\ is the density contrast between the two 

media). The cloud acceleration timescale, face, i-C- the amount of time it takes to accelerate the cloud to the velocity of the 
post-shock background Vps j = |V, is (e.g. Klein et al. 1994): 



_16pe^_16 

face — „ T7 ~ n VX^crush 

9 P i Vi 9 



100 kms- 



[yr]. 



(1) 
(2) 



The above equation show that it takes a long time to couple dynamically pre-shock dense H2 gas with the background flow 
of tenuous gas. The cloud crushing time (time for the inter-cloud shock to go through the cloud) to acceleration time ratio is 
femsh/facc = 9/16y'X' which shows that as soon as the density contrast is higher than w 10~, the cloud will be compressed on a 
timescale much shorter than the acceleration time. In the SQ shock, the high density contrast between cold clouds and the hot 
tenuous flow (x > 10^, see Guillard et al. 2009) would lead to deceleration timescales longer than the shock age. 
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On the other hand, in the model of H2 formation in the shock proposed by Guillard et al. (2009), the molecular gas is formed 
out of moderate density medium («h ~ 1 cm"^) that is compressed and decelerated while becoming molecular. Therefore, in this 
context, the coupling is much more efficient because the density contrast is much lower (x < 10^). This in-situ formation is a 
solution to explain the high velocity dispersion seen in the spectra, although we do not exclude the presence of pre-shock dense 
clouds. 



